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1 ^1 Abstract 

I In this paper we study the most-demanding predicate for computing the Euclidean Voronoi 

^ diagram of axes-aligned line segments, namely the Incircle predicate. Our contribution is two- 

^^ fold: firstly, we describe, in algorithmic terms, how to compute the Incircle predicate for axes- 

^^ aligned line segments, and secondly we compute its algebraic degree. Our primary aim is to 

^^ minimize the algebraic degree, while, at the same time, taking into account the amount of 

operations needed to compute our predicate of interest. 

r^ In our predicate analysis we show that the Incircle predicate can be answered by evaluating 

^•^ the signs of algebraic expressions of degree at most 6; this is half the algebraic degree we get 

I when we evaluate the Incircle predicate using the current state-of-the-art approach. In the most 

1^ demanding cases of our predicate evaluation, we reduce the problem of answering the Incircle 

. 1^ predicate to the problem of computing the sign of the value of a linear polynomial (in one 

Si^ variable), when evaluated at a known specific root of a quadratic polynomial (again in one 

}h variable). Another important aspect of our approach is that, from a geometric point of view, we 

answer the most difficult case of the predicate via implicitly performing point locations on an 

appropriately defined subdivision of the place induced by the Voronoi circle implicated in the 

Incircle predicate. 
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1 Introduction 

The Euclidean Voronoi diagrams of a set of line segments is one of the most well studied structures 
in computational geometry. There are numerous algorithms for its computation [6l [161 HBl [Ml 
[HI [H [T7]. These include worst-case optimal algorithms that use different algorithmic paradigms, 
such as the divide-and-conquer paradigm [24j or the sweep- line paradigm [8]. An interesting and 
efficient class of algorithms rely on the randomized incremental construction of the Voronoi diagram 
[H |T7]. From the implementation point of view, there are algorithms that assume that numerical 
computations are performed exactly |22| I14| . i.e., they follow the Exact Geometric Computation 
(EGC) paradigm [25], as well as algorithms that use floating-point arithmetic |12^ [23l 111) : the 
latter class of algorithms does not guarantee exactness, but rather topological correctness, meaning 
that the output of the algorithm has the correct topology of a Voronoi diagram. In terms of 
applications, these include computer graphics, pattern recognition, mesh generation, NC machining 
and geographical information systems (CIS) — see |16| [TSl [31 [Til E], and the references therein. 

Efficient and exact predicate evaluation in geometric algorithms is of vital importance. It has 
to be fast for the algorithm to be efficient. It has to be complete in the sense that it has to cover all 
degenerate cases, which, despite that fact that they are "degenerate" from the theoretical/ analysis 
point-of-view, they are commonplace in real world input. In the EGC paradigm context, exactness 
is the bare minimum that is required in order to guarantee the correctness of the algorithm. The 
efficiency of predicates is typically measured in terms of the algebraic degree of the expressions (in 
the input parameters) that are computed during the predicate evaluation, as well as the number 
(and possibly type) of arithmetic operations involved. The goal is not only to minimize the number 
of operations, but also to minimize the algebraic degree of the predicates, since it is the algebraic 
degree that determines the precision required for exact arithmetic. Degree-driven approaches for 
either the evaluation of predicates, or the design of the algorithm as a whole, has become an 
important question in algorithm/predicate design over the past few years [H [191 121 13 [13 13 I20j . 

In this paper we are interested in the most demanding predicate of the Euclidean Voronoi 
diagram of axes- aligned line segments. Axes- aligned line segments, or line segments forming a 45- 
degree angle with respect to the axes, are typical input instances in various applications, such as 
VLSI design |2H 110) . However, although the predicates for the Euclidean Voronoi diagram of line 
segments have already been studied [4j, the predicates for axes-aligned or ortho-45° input instances 
have not been studied in detail in the Euclidean setting. In the sections that follow, we analyze 
the Incircle predicate in this setting: given three sites Si, S2 and S3, such that the Voronoi circle 
y( Si, 5*2, 5*3) exists, and a query object O, we seek to determine if O intersects the disk D bounded 
by ^(^i, ^2, 5*3), touches D or is completely disjoint from D. In our context Si, S2, S3 and O are 
either points or axes-aligned (open) line segments. Our aim is to minimize the algebraic degree 
of the expressions involved in evaluating the Incircle predicate. We show that we can answer the 
Incircle predicate using polynomial expressions in the input quantities whose algebraic degree is at 
most 6. This is to be compared: (1) against the generic bound on the maximum algebraic degree 
needed to compute the Incircle predicate, when we impose no restriction on the geometry of the line 
segments, which is 40 |4J, and (2) against the specialization/simplification of the approach in [3], 
when we consider axes-aligned line segments. With respect to the latter case, our algebraic degrees 
are never worse, while in the most demanding case we have reduced the degree by a factor of two 
(see also Table II]). 

The rest of our paper is structured as follows. In Section [2] we give some definitions, compare 
our approach to that in [4j, and detail some of the tools that we use in the Incircle predicate analysis. 
In Sections [3][7] we describe how we evaluate the incircle predicate for different configurations of the 
sites Si, S2, S3 and O. In Section [8] we detail plans for future work. 
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Table 1: Maximum algebraic degrees for the eight types of the Incircle predicate according to: j^ 
for the general and the axes-aligned segments case, and this paper. Top/Bottom table: the query 
object is a point/segment. 



2 Definitions and preliminaries 

Given three sites 5i, 52, and 53 we denote their Voronoi circle by V{Si, 82,33) (if it exists). There 
are at most two Voronoi circles defined by the triplet (5i, 52, 53); the notation V{Si, S2, S3) refers 
to the Voronoi circle that "discovers" the sites Si, S2 and 53 in that (cyclic) order, when we walk 
on the circle's boundary in the counterclockwise sense. Given a fourth object O, which we call 
the query object, the Incircle predicate lncircle(5i, 52, 53, O) determines the relative position O with 
respect to the disk D bounded by V{Si, S2, S3). The predicate is positive if O does not intersect 
D, zero if O touches the boundary but not the interior of D, and negative of the intersection of O 
with the interior of D is non-empty. 

The Voronoi circle of three sites does not always exist. In this paper, however, we assume that 
the Incircle predicate is called during the execution of an incremental algorithm for computing the 
Euclidean Voronoi diagram of line segments, and thus the first three sites are always sites related to 
a Voronoi vertex in the diagram. Note that the value of the Incircle predicate does not change when 
we circularly rotate the first three arguments. In that respect, there are only four possible distinct 
configurations for the type of the Voronoi circle: PPP, PPS, PSS and SSS, where P stands for 
point and 5 stands for segment. For example, a Voronoi circle of PSS type goes through a point 
and is tangent to two segments. This gives eight possible configurations for the Incircle predicate, 
two per Voronoi circle type. 

The predicates for the Euclidean Voronoi diagram of line segments, in the context of an in- 
cremental construction of the diagram, have already been studied by Burnikel |3]. According to 
Burnikel's analysis the most demanding predicate is the Incircle predicate: assuming that the input 
is either rational points, or segments described by their endpoints as rational points, Burnikel shows 
that the Incircle predicate can be evaluated using polynomial expressions in the input quantities, 
whose algebraic degree is at most 40; this happens when the Voronoi circle is of SSS type and the 
query object is a segment (see also the line dubbed "General |5" in Table II]). Considering Burnikel's 
approach for the case of axes-aligned line segments, and performing the appropriate simplifications 
in his calculations, we arrive at a new set of algebraic degrees for the various configurations of the 
Incircle predicate (see line dubbed "Axes-aligned |3]" in TablefTl); now the most demanding case the 
is PPS case, which gives algebraic degree 8 and 12, when the query object is a point and a segment, 
respectively. 

In Sections |3]J7] we analyze, in more or less detail, all eight possible configurations for the Incircle 
predicate, and show how we can reduce the algebraic degrees for the PPS case from 8 and 12 to 



6. This is done by means of three key ingredients: (1) we formulate the Incircle predicate as an 
algebraic problem of the following form: we compute a linear polynomial L{x) = hx + ^o and a 
quadratic polynomial Q{x) = q2x'^ + qix + qo, such that the result of the Incircle predicate is the 
sign of L{x) evaluated at a specific root of Q{x), (2) for the PPS and PSS cases, we express the 
Incircle predicate as a difference of distances, instead of as a difference of squares of distances, and 
(3) we reduce the PPSP case to the PPPS case. Regarding the first ingredient, we describe in the 
following subsection how we can do better than finding the appropriate root of Q{x) and substitute 
it in L{x) (this is essentially what is done in |3]). Regarding the second and third ingredients we 
postpone the discussion until the corresponding sections. There is one final tool that we will be 
very useful in order to simplify our analysis: in order to reduce our case analysis we make extensive 
use of the reflection transformation through the line y = x; see Subsection |2.2| for the details. 

2.1 Evaluation of the sign of L(x) = hx+lo at a specific root oi Q(x) = q2x'^+qiX+qQ 

Let L(x) = Iix + Iq and Q{x) = q2x'^ +qix + qQ be a linear and a quadratic polynomial, respectively, 
such that Q{x) has non-negative discriminant. Let the algebraic degrees of /i, /2, 92, qi and go be 
5;, 5i + 1, (5g, 5q + 1, and 5q + 2, respectively. We are interested in the sign of -L(r), where r is one 
of the two roots xi < X2 of Q{x). In our analysis below we will assume, without loss of generality 
that li,q2 > 0. 

The obvious approach is to solve for r and substitute into the equation of L{x). Let Aq = 
qf — 4^290 be the discriminant of Q{x). Then r = {—qi =b y^ IS.q) / {2q2) , which in turn yields 
L{r) = {hqi + 2lQq2 ± y/~AQ) / {2q2) ■ Computing the sign of L{r) is dominated, with respect to the 
algebraic degree of the quantities involved, by the computation of the sign of liqi + 2/o'?2 i /i iVAg. 
Evaluating the sign of this quantity amounts to evaluating the sign of {hqi + 2/0^2)^ — ^i^Q, which 
is of algebraic degree 2(6i + dq + 1). 

Observe now that evaluating the sign of L(r) is equivalent to evaluating the sign of Q{x*), and 
possibly the sign of Q'{x*), where x* = —f- stands for the root of L{x). Indeed, if Q{x*) < 0, we 
immediately know that L{r) < if r = xi, or that L{r) > if r = X2. If Q{x*) > 0, we need to 
additionally evaluate the sign of Q'{x*) = 2q2X* + qi. If Q'{x*) < 0, we know that x* < xi,a:2, 
which implies that L{r) > 0, whereas if Q'{x*) > 0, we have x* > xi,X2, which gives L{r) < 0. 
Finally, if Q{x*) = 0, we still need to evaluate the sign of Q'{x*). If Q'{x*) < 0, then x* = xi, 
and thus L{r) = if r = xi, and L{r) > if r = X2- Similarly, if Q'{x*) > 0, then x* = X2, and 
thus L{r) < if r = xi, and L{r) = if r = X2- There is one last case to consider: Q'[x*) = 0. 
Given that Q{x*) = 0, this can happen only if xi = X2 = x*, in which case we deduce L(r) = 0. 
Since Q(x*) = {ifqo — WqiIq + 92^0)/^!' evaluating the sign of Q{x*) means evaluating the sign of an 
algebraic expression of degree 26i + 6q + 2. Moreover, Q'{x*) = {hqi — 2q2lo)/h; hence, evaluating 
the sign of Q'{x*) reduces to evaluating the signs of liqi — 2g2^o and /i, the degrees of which are 
^l + 5q + ^ and (5;, respectively. 

Notice that the latter among the two approaches described above is never worse than the first 
one; in fact, if 5g > it gives a lower maximum algebraic degree. We summarize this observation 
in the following lemma. 

Lemma 1. Let L{x) = Iix + Iq, li 7^ 0, andQ{x) = q2X^ + qix + qo, (72 7^ 0, 6e a linear and quadratic 
polynomial, respectively, such that the discriminant ofQ(x) is non-negative. If the algebraic degrees 
ofli, I2, Q2, Qi o.nd qQ be 5i, di + 1, Sq, Sq + 1, and 5q + 2, respectively, then we can evaluate the sign of 
L{r), where r is a specific root of Q{x), using expressions of maximum algebraic degree 25i + 5q + 2. 



2.2 Reflection transformation 



Let 7^ : E^ — )• E^ denote the reflection transformation tlirougli tlie line y = x. TZ maps a 
point (x,y) € E^ to the point {y.,x) G E^. The reflection transformation preserves circles and 
line segments and is inclusion preserving. This immediately implies that, given a Voronoi circle 
V{Si, S2, 5*3) defined by three sites Si, S2 and S3, and a query point Q, Q lies inside, on, or outside 
the Voronoi circle V{Si, S2, S^) if and only if 1Z{Q) lies inside, on, or outside the Voronoi circle 
V (TZ{S2) ,'R-{Si) ,1Z{S3)) (cf. Fig. [T]for the case where Si and ^2 are points and S^ is a segment). 
Hence, lncircle(S'i, 5*2, 53, Q) = lncircle(7^(52),7^(S'i), 7^.(53), 7^(Q)). Notice that reflection reverses 
the orientation of a circle, which is why we consider the Voronoi circle V {TZ{S2) ,'R'{Si) ,1Z{S3)) in- 
stead of the Voronoi circle V {lZ{Si) ,TZ{S2) ,TZ{Sz)) . The same principle applies in the case where the 
query object is a line segment QS: lncircle(S'i, 52, 5*3, QS) = \nc\rc\e{TZ{S2),'R'{Si),TZ{S3),lZ{QS)). 
As a final note, the reflection transformation IZ maps an x-axis parallel segment to a y-axis 
parallel segment, and vice versa. This property will be used, in the sections that follow, to reduce 
the analysis and computation of the Incircle predicate, where one of the SiS is y-axis parallel, to 
the case where one of the Si 's is x-axis parallel. 
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Figure 1: lncircle(A, B, CD, Q) is equivalent to lncircle(7^(5), 7^(^), 7^(CD), 7^(Q)), where 7^ stands 
for the image of / under the reflection transformation through the line y = x. 



3 The PPP case 

As of this section, we discuss and analyze the Incircle predicate for each of the four possible config- 
urations for the Voronoi circle. We start with the case where the Voronoi circle is defined by three 
points A, B and C. 



3.1 The query object is a point 

This is the well known Incircle predicate for four points A, B, C and Q, where Q is the query point, 
and it amounts to the computation of the sign of the determinant 
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Its algebraic degree is clearly 4. 



3.2 The query object is a segment 

Let QS be the query segment. In this case, we must first check that relative position of Q and S 
with respect to V{A,B,C) using lncircle(A, i?, C, /), / € {Q,5}. If at least one of Q and S lies 
inside V{A, B, C), we clearly have lncircle(A, B, C, QS) < 0. 

Otherwise, we have to examine if the segment QS intersects with V{A, B, C). This is equivalent 
to point-locating the points Q and S in the arrangement of the lines y = ymin, V = Umax and x = xk 



if QS is X-axis parallel or, x 



X and y = yx if QS is y-axis parallel, where x^. 



(resp., ymin ,ymax) are the extremal points of V{A, B, C) in the direction of the x-axis (resp., y-axis). 
In fact, the case where the segment QS is y-axis parallel can be reduced to the case where the query 
segment is x-axis parallel by noting that lncircle(^, B, C, QS) = \<nc\rc\e{Ti{B),'R{A),lZ{C),TZ{QS)) 



m 



(see Section 2.2). We will therefore restrict our analysis to the case where QS is x-axis parallel. 

We first determine if Q lies outside the band delimited by the lines y = ymin and y = ymax] 
this case we immediately get \nc\rc\e(A,B,C,QS) > 0. Otherwise, if Q lies inside the band (resp., 
Q lies on either y = ymin or y = ymax), we check the relative positions of Q and S against the line 
X = Xk] the segment QS intersects (resp., is tangent to) V{A,B, C) if and only if Q and S lie on 
different sides of the line x = xk- 

In order to determine the relative position of Q with respect to the lines y = ymin and y = ymax, 
we will evaluate a quadratic y-polynomial that vanishes at ymin and ymax'- let T{y) = t2y'^+tiy+to be 
this polynomial. Having computed this polynomial, yq £ {ymin, Umax) if and only if sign(T{yQ)) = 
-sign{t2), yq [ymin,y7nax] if and only if sign{T{yQ)) = sign{t2), and, finally, yg £ {ymin,ymax} 
if and only if sign{T(yQ)) = 0. 

To evaluate such a polynomial, we first observe that every point {x,y) on V{A,B,C) satisfies 
lncircle(A, B, C, (x, y)) = 0. Expanding the four-point Incircle determinant in terms of x, we end up 
with a quadratic polynomial U{x; y) = u^x^ -|- U\x -|- Mo(y) for Incircle(j4, S, C, (x, y)), where 
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For a fixed value y* of y, the roots of U{x] y*) are the points of intersection of the line y = y* with 
the Voronoi circle V{A, B, C). U{x; y*) has no real roots if y* [ymin, ymax], has two distinct roots 
if y* G (ymin, ymax) and has a double root if y* £ {ymin, ymax}- In the last case, the discriminant 
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4:U2Uo{y*) of U{x; y*) has to vanish. Now consider the discriminant as a polynomial 



of y. Clearly, A[/(y) is a quadratic y-polynomial, with a strictly negative, since the points A, B 
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Figure 2: Relative positions of the x-axis aligned query segment QS with respect to the lines x = xk, 
y ^ ymin, y ^ ymax- 
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and C are not collinear. Moreover, Aij{y) vanishes for y G {ymimymax}, hence it may serve as the 
quadratic polynomial T{y) we were aiming for. More specifically, T{y) := IS.u{y) = t2y'^ + tiy + to 
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In an analogous manner, we can evaluate a quadratic x-polynomial that vanishes at Xmin and 
Xmax, which we call S{x). More precisely, S{x) = S2X^ + s\x + sq, where S2 = — 4m|, s\ = —Au2Ui 
and So = ^1 + 4u2U3. In order to determine the relative position of Q and S with respect to the 
line X = xk, we use the fact that xk = ^{xmin + Xmax) = ~ f^- Hence, using the fact that S2 < 0, 
checking on which side of x = xk lies point /, for I G {Q,S}, amounts to determining the sign 
sign{xK — xj) = sign{2s2Xi + si). 

The algebraic degrees of uq, ui, U2, M3, and wi are 4, 3, 2, 3, and 3, respectively. Therefore, the 
algebraic degrees of ^2, ti, to, S2, si, and sq are 4, 5, 6, 4, 5, and 6, respectively. This implies that 
the algebraic degree of T[yQ) is 6, while the algebraic degree of S2X1 + S1, I G {Q, S}, is 5. We, thus, 
conclude that we can answer the Incircle predicate in the PPPS case by evaluating expressions of 
maximum algebraic degree 6. 

4 The SSS case 

In this section we consider the case where the Voronoi circle is defined by three axis-aligned seg- 
ments AB, CD and FG. In order for the circle V{AB,CD,FG) to be well defined, exactly two 
of these segments must parallel to each other, while the third perpendicular to the other two. 
Given that V{AB,CD,FG) = V{FG,AB,CD) = V{CD,FG,AB), we can assume without loss 
of generality that the first two segments are parallel to each other, and thus the third is per- 
pendicular to the first two. Hence we only have to consider two cases: (1) AB, CD are j;-axis 
parallel and FG is y-axis parallel, and (2) AB, CD are y-axis parallel and FG is x-axis parallel. 
In fact the second case can be reduced to the first one by noting that \r\c\rc\e{AB,CD,FG,Q) = 



\narc\e{TZ{CD),TZ{AB),n{FG),n{Q)) (see Section [Oj) . We shall, therefore, assume that AB, CD 
are x-axis parallel and FG is y-axis parallel. 

4.1 The query object is a point 

Let Q be the query point. Since the center K of V{AB,CD,FG) lies on the bisector of the lines 
Iab and icD, and the radius p of the circle is the distance of K from either Iab or icD (i-e., half 
the distance of the two lines), we have 

K = ix,+y^,y<^), p=\y^. (d 

To answer the Incircle predicate for Q, we first examine if Q and K lie on the same side with 
respect to the lines iAB, f-CD and ipG- If this is not the case, we immediately conclude that 
\r\c\rc\e{ AB,CD,FG,Q) > 0. Otherwise we must compare the distance d{Q,K) of Q from K 
against the Voronoi radius p. More precisely: \'nc\'cc\e{AB,CD,FG,Q) = sign{(P{Q,K) — p^), 
where A{(f{Q,K) - p^) = 4(xf - xq){1 + yc - Va) + {vc + VA - 2yQ)^, which is an algebraic 
expression of degree 2 in the input quantities. Given that the sideness tests for Q against the lines 
^AB^ ^CD and ipc are of degree 1, we conclude that answering the Incircle predicate in the SSSP 
case amounts to computing the signs of expressions of algebraic degree at most 2. 

7 



4.2 The query object is a segment 

Let QS be the query segment. We first determine if tlie endpoints Q and/or S of QS lie inside 
V{AB, CD, FG), in which case we immediately get lncircle(Ai?, CD, FG, QS) < 0. Otherwise, we 
must consider the orientation of QS and make the appropriate checks. 

Assume first that QS is x-axis parallel. We first check if Q is inside the band By delimited by 
the lines Iab and IcD- If Q lies outside By, we immediately get that lncircle(^i?, CD, FG, QS) > 0. 
Otherwise, we have to determine the relative positions of Q and S with respect to the line x = xk, 
where xk = xp + ^(yc' ~ Va), by evaluating the signs of xq — xk and xs — xk- If Q lies inside 
By (resp., on the boundary of By, QS intersects (resp., is tangent to) V{AB, CD, FG), if and only 
if Q and S lie on different sides of the line x = xk, i-e., if and only if (xq — xk){xs — xk) < 0. 
Determining if Q lies inside By amounts to computing the signs of yg — yA and yq — yc, which are 
degree 1 quantities. The quantities xq — xk and xg — xk are also of degree 1, which implies that 
we can answer the Incircle predicate in this case using quantities of algebraic degree up to 2 (the 
algebraic degree needed to evaluate lncircle(^i?, CD, FG,/), / G {Q-,S} dominates the degrees of 
all other quantities to be evaluated). 

Consider now the case where QS is y-axis parallel. We first need to check if the line (.qs, 
intersects with V{AB, CD, FG). To do this we need to evaluate the sign of quantity \xq — xk\ — P, 
where p is given by (IT| . Computing the signs of xq— xk and yc — yA, we may express | xq — x;^ | — /> as 
a polynomial expression in the input quantities; its algebraic degree is, clearly, 1. If \xq — xk\—P > 0, 
£qs does not intersect V{AB,CD,FG), and we immediately get \nc\rc\e{ AB,CD,FG,QS) > 0. 
Otherwise, if \xq — xk\ — p < (resp., \xq — xk\ — p = 0) iqs either intersects with (resp. either is 
tangent to) the Voronoi circle or does not intersect the Voronoi circle at all. To distinguish between 
these two cases we have to determine if the points Q and 5 lie on different sides of the line y = yx- 
the segment QS intersects with (resp., is tangent to) the Voronoi circle V{AB, CD, FG) if and only 
if {yQ-yK){ys-yK) < O. Since yK = l{yA + yc) (see rel. ([I|), determining the signs sign{yQ-yK) 
and sign{ys — yx) amounts to computing the sign of quantities of algebraic degree 1. As in the case 
where QS is x-axis parallel, the algebraic degree for evaluating the Incircle predicate is dominated 
by the algebraic degree for evaluating lncircle(y4B, CD, FG, I), I £ {Q, S}, which is 2. 

5 The generic approach for the evaluation of the Incircle predicate 
in the PPS and PSS cases 

In this section we present our approach for evaluating the Incircle predicate in a generic manner. 
The approach presented is applicable when the Voronoi circle is defined by at least one point and 
at least one segment, i.e., we can treat the cases PPS and PSS. 

Let K = {xK^yx) be the center of the Voronoi circle defined by the sites Si, 52, 5*3, that touches 
the sites ^i, S2 and S^ in that order when we traverse the Voronoi circle in the counterclockwise 
sense. As already stated, we want to evaluate the Incircle predicate for a query point or a query 
line segment with respect to this circle. To do this we compute a quadratic polynomial P{x) that 
vanishes at xk^ while using geometric considerations and the requirement on the orientation of the 
Voronoi circle, we can determine which of the roots xi < X2 of P{x) corresponds to xk- Regarding 
yK, the situation is entirely symmetric. We also compute a quadratic polynomial T{y) that vanishes 
at yK and, as for xk, we can determine which of the two roots yi < 2/2 of T(y) corresponds to yK- 
Moreover, in all cases xk and yK are linearly dependent, which means that we may express yK as 
yK = ^xk + ^, where ai, oq and /3 are polynomials in the input quantities. 



5.1 The query site is a point 

Let Q be the query point. Since at least one of 5*1, 5*2 and S^ is a point A, determining the Incircle 
predicate amounts to evaluating the sign of the quantity (P{K, Q) — d?{K, A) = {xk — xq)"^ + {yK — 
yQf'~{xK — XAf' — {yK—yAf'- Replacing yi^, using the relation yx = ^^K + ^, and gathering the 
terms of xk, we get lncircle(S'i, S2, S3, Q) = ^{hxR + h), where h = 2I3{xq - xa) + 2ai(yQ - va) 
and Jo = I3{xq + y'n — x\ — y\) — 2ao{yQ — yA)- If h = 0, the we can immediately evaluate the 
Incircle predicate by evaluating the signs of Iq and /3. Otherwise, deciding the Incircle predicate 
reduces to evaluating the sign of /?, as well as the sign of Iix + Iq, evaluated at a specific known 
root of a quadratic polynomial P{x) = p2x'^ + pix + pq (it is the root of P{x) that corresponds to 



Xk)- This is exactly the problem we analyzed in Subsection 2.1 

Let us now analyze the algebraic degrees of the expressions above. As we will see in the upcoming 
sections (see Sections [6] and l7| , P{x) is a homogeneous polynomial in terms of its algebraic degree. 
Letting 6x the algebraic degree oi p2, the algebraic degrees of pi and pQ become Sx + 1 and 6x + 2. 
Let also 6a be the algebraic degree of ai. In our context, the algebraic degree of oq is always one 
more that the degree of ai, i.e., it is 6a + 1, whereas the algebraic degree of /3 is always equal to 
that of ai. This implies that the algebraic degrees of Ii and Iq are (5q, + 1 and 6a + 2, respectively. 
Applying Lemma [T] we conclude that we can resolve resolve the Incircle predicate using expressions 
of maximum algebraic degree 2{6a + 1) + 6x + 2 = 26a + 6x + 4:. 

5.2 The query site is a segment 

Let QS be the query segment. The first step is to compute lncircle(S'i, 52, 5*3, Q) and, if needed, 
lncircle(5i, 52, S's, 5). If at least one Q and S lies inside the Voronoi circle V{Si, S2, S3), we get 
lncircle(5i, S2, S3, QS) < 0. Otherwise, we need to determine if the line iqs intersects V{Si, S2, S3). 
If Iqs does not intersect the Voronoi circle, we have lncircle(5'i, 5*2, S3, QS) > 0. If igs intersects 
the Voronoi circle we have to check if Q and S lie on the same or opposite sides of the line ^qs{K) 
that goes through the Voronoi center K and is perpendicular to Iqs- Notice that since QS is axes- 
aligned, the line £Qg{K) is either the line x = xk or the line y = yx- Since at least one of Si, S2 
and 5*3 is a segment CD, answering the Incircle predicate is equivalent to comparing the distance of 
K from the line iqs to the segment CD: 

lncircle(5i,52,53,£Qs) = d{K,iQs) - d{K,CD). (2) 

We can assume without loss of generality that CD is a;-axis parallel, since, otherwise we can re- 



duce lncircle(S'i, 52,53,(35') to lncircle(7^(52),7^(5l),7^(53), TZ{QS)) (see Section [2^, in which case 
1Z{CD) is X-axis parallel. Let us now examine and analyze the right-hand side difference of ^. 

Assume first that the segment QS is rc-axis parallel. In this case the equation of i.Qs is y = yq, 
and, hence, d{K, iqs) = \yK ~yQ\- Recall that yx is a specific root of a quadratic polynomial T{y). 
Therefore, determining the sign of yx — yq reduces to evaluating the sign of T{yq) and T'{yq). Let 
T(y) = ^22/^ + tiy + to be this polynomial, and let 6y, 6y + l, 6y + 2 be the algebraic degrees of t2, ti 
and toi respectively (as for P{x), T[y) is a homogeneous polynomial). Consider now the case where 
QS is y-axis parallel. The equation of iqs is 3; = xq, and, hence, d{K, iqs) = \xk — xq\. As in the 
X-axis parallel case, xk is a specific known root of the quadratic polynomial P{x), and determining 
the sign of xk — xq amounts to evaluating the sign of P{xq) and P'{xq). Last but not least, since 
the segment CD is x-axis parallel, d{K, CD) = {yx — yc\- As before, we can determine the sign of 
yK — yc by evaluating the signs of T{yc) and T'{yc)- 

Having made the above observations, we conclude that, if QS is x-axis parallel, 

lncircle(5i,52,53,£QS') = \yK - yq\ - \yK - yc\ = Jwk + Jo, 



where Ji and Jq are given in the following table. 



VK -VQ 


VK-yc 


Ji 


Jo 


>o 


>o 





yc-yq 


<o 


2 


-VQ - yc 


<0 


>o 


-2 


yq + yc 


<o 





-yc + yq 



Clearly, if Ji = we have lncircle((S'i, 52, S'3,£q5) = sign{Jo). Otherwise, given that yx is a root 
of T{y), evaluating lncircle(5i, 52, S'3,£q5) can be done using the analysis in Subsection 2.1 Since 



the algebraic degrees of Ji and Jq are and 1, respectively, we deduce, by Lemma [T] that we can 
resolve the Incircle predicate using expressions of algebraic degree at most 2 ■ + 5y + 2 = 6y -\- 2. 

For the case where QS is y-parallel we use the fact that yx = ^xk + tt- Using this linear 
dependence between xk and y^j we get 

lncircle(5i,52,53,£Q5) = \xk - xq\ - \yK - yc\ = 

where Li and Lq are given in the following table. 



/3 



{LiXK + Lq), 



XK - XQ 



> 



<0 



yx -yc 



Li 



>0 



-ai + P 



<0 



ai +/3 



> 



-ai — j3 



Ln 



/3(yc - xq) - ao 



P{-yc -xq) + ao 



P{yc + xq) - ao 



P{-yc + XQ) +ao 

Otherwise, given that xk is a known root 

As in the previous 



< I Qi - /3 

If Li = 0, \ndrc\e{Si,S2,S3,eQs) = sign{Lo)sign{l3). 

of P{x), determining the sign of Lixk + Lq can be done as in Subsection 2.1 

subsection, we let 6a be the algebraic degree of ai (and also of f3), which means that the degree of 

ao is 6a + 1. Hence, the algebraic degree of Li is 6a, whereas that of Lq is maxj^Q, + 1,1} = 6a + 1- 

By Lemma [T] in order to evaluate the sign LiXk + Lq we need to compute the signs of expressions 

of algebraic degree at most 26a + 6x + 2. 

As we mentioned at the beginning of this subsection, if lncircle(5i, ^2, S^jIqs) < 0, we need to 
check the position of Q and 5 with respect to the either line x = xk (if QS is x-axis parallel), 
or the line y = yx (if QS is y-axis parallel). To check the position of /, / G {Q,S}, against the 
line X = xk, we simply have to compute the signs of P{xj) and P'[xi). The algebraic degrees of 
these quantities are 6x + 2 and 6x + 1, respectively. In a symmetric manner, to check the position 
of /, / G {Q,S}, against the line y = yx, we simply have to compute the signs of T{yj) and 
T'{yj). The algebraic degrees of these quantities are 6y + 2 and 6y + 1, respectively. Notice that 
in both cases for the orientation of QS, the algebraic degree of the quantities whose sign needs 
to be evaluated to resolve the Incircle predicate are never greater than those computed above for 
evaluating lncircle(Si, ^2, 53,^q5). Recalling that, in order to evaluate lncircle(5i, £'2; 'S'3,<35'), the 
first step is to evaluate lncircle(S'i, ^2, S3, Q), and, if needed, lncircle(S'i, ^2, 5*3, S), we conclude that 
in order to evaluate the Incircle predicate when the query object is a segment we need to compute 
the sign of polynomial expressions of algebraic degree at most max{26a + 6x + 4:,6y + 2}. 



6 The PPS case 

Let A and B be the two points and CD be the segment defining the Voronoi circle. With- 
out loss of generality we may assume that CD is x-axis parallel, since otherwise we can reduce 



\ncnc\e{A,B,CD,Q) to \nchc\e{n{B),TZ{A),TZ{CD),n{Q)), as described in Section 2.2 
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degrees of p2, Pi and po are 1, 2 and 3, respectively. Moreover, in this case i/k 
where the algebraic degrees of ai, ao and /3 are 1, 2 and 1, respectively (i.e., 6a ■ 



6.1 The query object is a point 

Let Q be the query point, and K be the center of V{A, B, CD). As we will see in the next subsection, 
the x-coordinate of iT is a root of a quadratic equation P{x) = p2X^ +pix + poi where the algebraic 

Sx = 1). By 

Subsection 5.1 we can evaluate lncircle(^, B, CD, Q) using algebraic expressions of maximum degree 
2-1 + 1 + 4 = 7. Below, we are going to show how to lower this maximum algebraic degree to 6. 

Clearly, for the Voronoi circle V{A, B, CD) to be defined, both A and B must be on the same 
side with respect to icD- Consider now Q: if Q does not lie on the side of icD that A and B lie, we 
have lncircle(^, B, CD, Q) > 0. Testing the sideness of /, I G {A, B, Q}, against icD simply means 
testing the sign of yj — yc, which is a quantity of algebraic degree 1. 

Suppose now that Q lies on the same side of icD as A and B, and assume, without loss of 
generality, that Orientation(A, C, -D) > (the argument in the case Orientation(A, C, L>) < 0, or 
when one of A and B lies on icD , is analogous) . Consider the result a of the orientation predicate 
Orientation(74, B, Q). In the special case o" = (i.e., Q lies on the line Iab), we observe that Q lies 
inside the Voronoi circle V{A, B, CD) if and only if Q lies on Iab and between A and B. This can 
be determined by evaluating the signs of differences xq — xa and xq — xb , which are both quantities 
of algebraic degree 1. 

If o" 7^ 0, we are going to reduce lncircle(A, B, CD, Q) to lncircle(A, B, Q, CD) (see also Fig. [3]). 






D 



Figure 3: Reducing \nc\rc\e{A, B,CD,Q) to \nc\rc\e{A, B , Q , C D) . Top/Bottom row: Q lies to the 
left/right of the oriented line Iab- Left/Right column: Q lies inside/outside V{A,B,CD). The 
dotted circle is the Voronoi circle of A, B and Q. 
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Suppose first that o" > 0, i.e., Q lies to the left of the oriented line Iab- Since A, B and CD appear 
on V{A, B, CD) in that order when we traverse it in the counterclockwise sense, we conclude that Q 
lies inside V{A, B, CD) (resp., lies on V{A, B, CD)) if and only if the circle defined by A, B and Q, 
does not intersect with (resp., touches) the line icD- To see this, simply "push" the Voronoi circle 
towards Q, while keeping its center on the bisector of A and B. Hence, \nc\rc\e{A,B,CD,Q) = 
— Incircle(^, -B, Q, CD). In a similar manner, if o" < 0, i.e., Q lies to the right of the oriented line 
^AB, Q hes inside V{A, B, CD) (resp., lies on V{A, B, CD)) if and only if the circle defined by A, 
B and Q intersects the line icD- Hence, \nc\rc\e{A, B , C D , Q) = \nc\rc\e{B,A,Q,CD). 

Summarizing our analysis above, we first need to determine on which side of icD Q hes: this 
a degree 1 predicate. If needed, the next step is to compute Orientation (yl, i3, Q), which is a 
degree 2 predicate. If Oner\tat\on(A, B , Q) = we need two additional tests of degree 1 to answer 
\nc\rc\e{A,B,CD,Q); otherwise, we observe that 



\nc\rc\e{A,B,CD,Q) 



-lncircle(yl, B, Q, CD), if Orientation(^, B,Q)>0 
lncircle(5, A, Q, CD), if Orientation(^, B,Q) <0 



As per Section 3.2 \nc\rc\e{A, B , Q , C D) or \nc\rc\e{B , A, Q , C D) can be answered using quantities 
of algebraic degree at most 6. 



6.2 The query object is a segment 



For this case we are going to follow the generic analysis presented in Section 5.2 Let QS be 



the query segment, and let K be the center of V{A,B,CD). K is an intersection point of the 
bisector of A and B and the parabola with focal point A and directrix the supporting line (.cd 
of CD. Solving the corresponding system of equations we deduce that, in the general case where 
A and B are not equidistant from Iqd (i-e., if Va 7^ yB)^ the x-coordinate of the Voronoi center 
xk: is a root of the quadratic polynomial P{x) = p2x'^ + pix + po, where p2 

iVB - yc){xA - Xb) - 2XBP2, Po = P2X% + {vc 



VB 



Pi 



'JC 



VB, 



[(4 



Xa) + ivA - yc)p2\ 



- yA / 0, 

while the 



y-coordinate of the Voronoi center hk, is a root of the quadratic polynomial T{y) = t2y + tiy + to, 



where t2 

XB 

VK 



H2yl + 2y% 



4(ys - yA), 

^yl + 



h = 

{xA 



^2yc 
^2^ 



yA 



xbd + {yi 



yB){xB - xa)'^ + 4:{yB - yA){y\ - y|), to = {xa 
-y%?- 



Moreover, yx and xk are linearly dependent: 



The 



^XK + ^, where ai = 2{xa - xb), a^ = x% + y%- x\- y\ and ^ = 2{yB - yA] 
roots xi < X2 of the polynomial P{x) (resp. yi < 2/2 of T{y)) correspond to the centers of the two 
possible Voronoi circles V{A,B,CD) and V{B,A,CD). The roots of P{x) or of T{y)) of interest 
are shown in the following two tables. 



Relative positions of A, B and CD 


Root of P{x) of interest 


yc <yA< ys 


Xl 


yc <yB < VA 


X2 


ys <yA< yc 


X2 


yA <yB < yc 


Xl 



Relative positions of A, B 


Root of T{y) of interest 


XA < XB 


2/2 


XA > XB 


yi 



The degrees of p2, Pi, Po, t2, ti and to are 1, 2, 3, 2, 3 and 4, respectively. Furthermore, the 
degrees of ai, ao and j3 are 1, 2 and 1, respectively. Applying the analysis in Subsection 5.2 (where 
^a = ^x = li Sy = 2), we deduce that we can answer the Incircle predicate using expressions of 
algebraic maximum algebraic degree max{2 • 1 + 1 + 2, 2 + 2} = 5. 
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Ki' 



"* A 

/ \ 
I 



•K2 



/ 




D 



-A. — -A. J JC — -X-j 

Figure 4: Voronoi circle defined by the point A and the hne segments CD and FG. Left: C-D, -FG 
are x-axis parallel. Right: CD is x-axis parallel and FG is y-axis parallel. 



For the special case ha 
^2 



yB, we easily get xk = ^{xa + xb) and Vk = ^, where U2 = 
{xb — xaY + ^{Va ~ Vc)^ ^1 ~ '^{VA — Vc)- In this case, if QS is x-axis parallel, we need to 
determine the sign of the quantity d(K,iQs) — d{K, CD) = |yx — UqI — IVK — 2/c|i or, equivalently, 
the sign of the quantity \U2 — Uii/qI — \U2 — Uiyc\, which is of algebraic degree 2. If QS is y-axis 
parallel, we need to evaluate the sign of the quantity d{K, iqs) — d{K, CD) = \xk — xq\ — {yx —yc\i 
or, equivalently, the sign of the quantity | f/i (x^ -|-xb — 2xq) | — 2| C/2 — C/iyg | , which is also of algebraic 
degree 2. Given, that the algebraic degree for the PPSP case is 6 (see previous subsection), we 
conclude that we can answer the Incircle predicate in the PPSS case by computing the signs of 
expressions of algebraic degree at most 6. 



7 The PSS case 



7.1 The query object is a point 

In this section we consider the case where the Voronoi circle is defined by two segments, a point 
and the query object is a point. Let A, CD and FG be the point and the two segments defining the 
Voronoi circle and let Q be the query point. Since each of CD, FG may be x-axis or y-axis parallel 
we have four cases to consider: (1) CD and FG are x-axis parallel, (2) CD and FG are y-axis 
parallel, (3) CD is x-axis parallel and FG is y-axis parallel, and (4) CD is y-axis parallel and FG 
is X-axis parallel. However, Cases (2) and (4) reduce to Cases (1) and (4), respectively, by simply 
performing a reflection transformation through the line y = x (see Section [2.2^ . More precisely, 
in both cases we have \nc\rc\e{A,CD,FG,Q) = \nc\rc\e{n{A),n{FG),n{CD),n{Q)). Thus, for 
Case (2), 1Z{CD) and TZ{FG) are x-axis parallel, while, for Case (4), 1Z{CD) is x-axis parallel and 
1Z{FG) is y-axis parallel. Therefore it suffices to consider Cases (1) and (3). In what follows, we 



follow the generic procedure described in Subsection 5.1 and refer to the notation introduced there. 



7.1.1 CD and FG are x-axis parallel 

We first notice that if Q does not lie inside the band B^ delimited by the IcD and ipGi it cannot 
be inside the Voronoi circle V{A,CD,FG). This can be easily checked by evaluating the signs of 
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VQ ~ yc aiid hq — i/F, which are quantities of algebraic degree 1. Suppose now that Q is inside B^ 
and notice that A has to he in Bx in order for the Voronoi circle V{A, CD, FG) to exist. 

Let K be the center of V{A,CD,FG). The y-coordinate of iC is, trivially, yx = 2^110 + yF), 
whereas the radius p of the Voronoi circle is equal to p = ^Ivc ~ VfI- Given that A is a point on 
V{A, CD, FC), we have that (i^{K, A) = p^ ■ Using the expressions for yx and p, we deduce that xk 
is a root of the polynomial P{x) = x^ +pix+pQ, where pi = 2xa and po = x\ + {yA — yc){yA — yF)- 
If xi < X2 are the two roots of P{x), the root that corresponds to xk is given in the table below 
(see also Fig. gleft)). 



Relative positions of A and CD 


Root of P{x) of interest 


yA > yc 


X2 


yA <yc 


Xl 



Moreover, in this case we have ai = 0, ao = 2/c + 2/F and /3 = 2. Therefore, the algebraic degrees 



involved in the evaluation of the Incircle predicate are 6a = ^x = 0. As per Subsection 5.1 the 
\nc\rc\e{A,CD,FC,Q) predicate can be evaluated using algebraic expressions of maximum degree 
2-0 + + 4 = 4. 

7.1.2 CD is X-axis parallel and FC is y-axis parallel 

The lines icD and ipc subdivide the plane into four quadrants Ri, R2, R3 and i?4. The bisector 
of i?i and i?3 is the line £1^3 with equation y = x + yc — xp, whereas the bisector of R2 and i?4 is 
the line £2,4 with equation y = —x + yc + xp. 

The center K of the Voronoi circle V{A, CD, FC) lies on both the bisector of icD and ipG^ as 
well as on the parabola that is at equal distance from A and icD] the equation of the latter is: 

{x - xa)^ - {yA - yc){'^y - yA - yc) = 0- (3) 

Assuming that A lies in Ri U i?3, the bisector of icD and ipc is ^1.3 ■ Substituting y in terms of 
X, using the equation of ^1^3, we deduce that the x-coordinate xx of i^ is a root of the quadratic 
polynomial P{x) = x'^+pix+po, where pi = 2{yc-yA-XA), andpo = {yc-yA)'^+XA-2xF{yc-yA)- 
Similarly, if A lies in R2U R4, xk is a root of the quadratic polynomial P{x) = x^ +pix +poi where 
pi = 2{yA -yc - xa), Po = {yc - yA)'^ + Xa + '^XF{yc - yA)- If xi < X2 are the two roots of P{x), 
the root that corresponds to xk is the same as in the case where FC is x-axis parallel. Moreover, 
in this case we have ai = 1, ao = yc — xp, (3 = 1, if A £ RiU R3, and ai = —1, ao = yc + xp, 
/3 = 1, if A £ R2L) R4. In both cases, the algebraic degrees involved in the evaluation of the Incircle 



predicate are 5a = Sx = 0. Again, as per Subsection 5.1 the \nc\\rc\e{A,CD,FG,Q) predicate can 



be evaluated using algebraic expressions of maximum degree 2-0 + + 4 = 4. 

7.2 The query object is a segment 

Let QS be the query segment, while the Voronoi circle is defined by the point A and the segments 
CD and FC. Let K = {xk, yx) be the center of the Voronoi circle. As in the previous subsection, 
it suffices to consider the cases where, either both CD and FC are x-axis parallel, or CD is x-axis 
parallel and FC is y-axis parallel. Recall that, in both cases, we have shown that xk is always a 
root of a quadratic polynomial P{x) = x^ +pix +poi where the algebraic degrees of pi and po are 
1 and 2, respectively. 
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7.2.1 CD and FG are x-axis parallel 

If QS is also X-axis parallel we first need to determine if QS lies inside the band B^ delimited by ^CD 
and ^FG- This is easily done by checking if Q lies inside B^, which in turn means checking the signs of 
yq—yc and yq — yp, as described in the previous subsection. Clearly, if Q is not inside the band B^, 
then lncircle(74, CD, FG, QS) > 0. Assume now that QS lies inside B^- The first step is to evaluate 
the lncircle(^,CZ),FG,Q) and, if necessary, \nc\rc\e{A,CD,FG,S). If lncircle(^, CL>, FG, Q) < 
or \nc\rc\e{A,CD,FG,S) < 0, then we immediately know that \ncnc\e{ A, C D, FG,QS) < 0. 
Otherwise, we simply need to determine on which side of the line x = xk Q and S lie: QS 
intersects the Voronoi circle V{A, CD, FG) if and only if Q and S lie on different sides of x = xk- 
Determining the side of x = xk on which the point I, I G {Q,S}, lies is equivalent to computing 
the sign of the difference xk — xj. This, in turn, reduces to computing the signs of the expressions 
P{xi) and P'{xj), which are expressions of algebraic degree 2 and 1, respectively. 

In the case where QS is y-axis parallel, we proceed according to the generic approach presented 
In this case yx = ^{yc + I/f), i-e., ai = 0, uq = yc + yp and j3 = 2. Moreover, 



5.2 



in Subsection 

T[y) is a linear polynomial T{y) = 2y — [yc + I/f)) thus the algebraic degrees of T(yj) and T'{yi) 
I G {Q, S"}, are 5y + l = 1 and 5y = 0, respectively. By applying the analysis of Subsection |5.2[ with 
^x = ^a = ^y = 0) we conclude that we can answer the Incircle predicate by evaluating the signs of 
expressions of algebraic degree at most max{2 • + + 4, + 1} =4. 

7.2.2 CD is X-axis parallel and FG is y-axis parallel 



For the purposes of resolving this case, we are going to follow the analysis of Subsection 5.2 
In the previous subsection we argued that in this case the center K = {x^^yx) of the Voronoi 
circle V{A, CD, FG) lies on the intersection of the parabola with equation (pi) and either the line 
y = X + yc — xp (if A G i?i U R^) or the line y = —x + yc + xp (if ^ S i?2 U R4). Solving in 
terms of y we deduce that yx is a root of the quadratic polynomial T{y) = y^ + tiy + tg, where 
ti = -2{yA + XA + XP - 2yc), to = [xA + xp)'^ + Va ~ '^VcixA + xp), if A G i?i U R^, whereas 
ti = 2{xA — yA — Xp), to = {xA — xp)"^ + yA~ "^XAyc + '^ycxp, if A G R2U R^. Notice that in both 
cases the algebraic degrees of ti and to are 1 and 2, respectively. Furthermore, if yi < y2 are the 
two roots of T{y), the root of T{y) of interest is given in the following table (see also Fig. Efright)). 



Relative positions of A and FG 


Root of r(y) of interest 


Xa > Xp 


2/2 


Xa < Xp 


yi 



Finally, as already described in the previous subsection, in this case we have ai = 1, uq = yc — xp, 
13 = 1, if A £ Ri U R3, and ai = —1, ao = yc + xp, (3 = 1, if A £ R2 U R4. We are now 



ready to apply the analysis of Subsection 5.2 with 6a = Sx = 6y = 0. We thus conclude that 
the predicate lncircle(A, CD, FG, QS) can be evaluated using algebraic quantities of degree at most 
max{2-0 + + 4,0 + 2} = 4. 

8 Conclusion and future work 

In this paper we have studied the Incircle predicate involved in the computation of the Euclidean 
Voronoi diagram for axes-aligned line segments. We have described in detail, and in a self-contained 
manner, how to evaluate this predicate. We have shown that we can always resolve it using poly- 
nomial expressions in the input quantities that are of maximum algebraic degree 6. 
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Our analysis is thus far theoretical. We would like to implement the approach presented in this 
paper and compare it against the generic implementation in CGAL [13]. Finally, we would like 
to study the rest of the predicates involved in the computation of the Voronoi diagram, as well as 
consider the ortho-45° case, i.e., the case where the segments are allowed to lie on lines parallel to 
the lines y = x and y = —x. 
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